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Abstract: C-parity-even quarkonia 7/(, jC and Xb,c with spin and 2 are produced via 
r~| ! two-gluon fusion. The expected cross section of the inclusive production of the quarkonia 

at the LHC, times the branching ratios of convenient decays, is up to tens of nanobarn 
per unit rapidity in the case of charmonia and around one nanobarn for the bottomonia. 
Measuring the quarkonia production as function of rapidity will allow to determine the 
gluon distribution function in nucleons in a very broad range of the Bjorken x from x ~ 
1CP 2 where it is already known, down to x ~ 10 -6 where it is totally unknown. The 
scale of the gluon distribution found from such measurements turns out to be rather low, 
Q 2 ~ 2.5 - 3GeV 2 , for char monia and rather large, Q 2 ~ 20GeV 2 , for bottomonia. We 
evaluate the scale by studying the next-to- leading-order production cross sections. 
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1 Introduction 

The small- x gluon distribution function in nucleons at a relatively low momentum scale 
is a fundamental quantity in high energy physics, determining the bulk of the collision 
processes. Apart from being of practical importance for evaluating the rate of many pro- 
cesses at high energies and of the background for new physics, the gluon distribution in 
nucleons has its own fundamental value as it collects many fine and subtle features of 
Quantum Chromodynamics. At a relatively low momentum scale and small x one expects 
the transition from the hard DGLAP regime to the soft nonperturbative pomeron [1] but 
their interplay is not fully understood. Theoretical models predict various gluon distri- 
bution functions g(x,Q 2 ), therefore, knowing it one can discriminate between the models. 
However, the experimental knowledge of this fundamental quantity is so far limited. 

At present, the low-x global parton analysis is based mainly on the deep inelastic scat- 
tering HERA data where quark (and antiquark) but not gluon distributions are measured 
directly. The gluon parton distribution functions (PDFs) are extracted from the derivative 
dF^ix, Q 2 ) /din Q 2 using the DGLAP evolution equation. For this reason the accuracy in 
the determination of the gluon densities is not too good. Moreover, in the range of very 
small x < 10~ 3 and at low momenta scale Q 2 ~ 2 — 3 GeV 2 the present-day gluon distri- 
butions are actually given by ad hoc extrapolations from the larger x data since this range 
has not been accessible by the previous data. 

In Fig. 1 we plot the low-x extrapolations from the CT10 [2] and NNPDF [3, 4] gluon 
distributions. The drop of the gluon flux x g (x, 2.5 GeV 2 ) at very small x is counter- 
intuitive: on the contrary, one expects that it should be roughly a constant, which would 
correspond to a constant cross section for minijet production, or even rise as a small 
power of 1/x, see Fig. 1, left. The unexpected behaviour of g(x) may be a result of 
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neglecting power and absorptive corrections that are probably non-negligible at relatively 
low Q 2 ~ 2.5 GeV 2 . It should be noted that the MSTW low-x NLO gluon distribution [5] 
becomes even negative at this low scale, which gives the idea of the uncertainty in the 
present-day knowledge. 



xg(x, 2.5 GeV 




xg(x, 20GeV 2 ) 
20 



log(x) 




log(jT) 



Figure 1. Gluon distribution function (times x) x g(x, Q 2 ) for the scales Q 2 = 2.5 GeV 2 appropriate 
for charmonia production (left), and Q 2 = 20 GeV 2 appropriate for bottomonia production (right). 
The shadowed areas are spanning the extrapolations of the CT10 (upper side) and of the NNPDF 
(lower side) NLO parton distributions. The curves show our extrapolations to the small- x range 
assuming xg(x) ~ const, (dot-dashed lines), xg(x) ~ l/x 0A (solid lines) and xg(x) ~ 1/x ' 2 
(dashed lines). On the right, the dashed line shows the extrapolation xg(x) ~ 1/cc ' 24 . The plots 
give the idea of the vast uncertainty in the present-day knowledge of the gluon distribution at very 
small x. 

The much higher energy of the LHC and a relatively low mass of the rj c and Xc mesons 
allows to probe the gluon distribution directly down to a few units of 10~ 6 . Indeed, the 
rjc with spin and Xc mesons with spin J = 0, 2 having positive C-parity are produced in 
the leading order (LO) via the simple gluon-gluon fusion gg — > r] c ,XcO,Xc2, and similarly 
for the bottomonia. The two-gluon fusion into spin-1 mesons such as J and Xci is 
forbidden by the Landau- Yang selection rule [6], therefore the rj c (0 h , 2980), Xco(0 ++ , 3415) 
and Xc2(2 ++ , 3556) mesons are, in this sense, privileged. A survey of quarkonia production 
in high energy collisions can be found in Ref. [7]. 

In the LO the inclusive production cross section of C-even quarkonia, integrated over 
the transverse momentum of a meson is given by a simple factorized equation [8] 1 

da(pp —7- quarkonium) 



dY 



xig(x\, hf) X2 g(x2, I^f) &(gg — > quarkonium) 



(1.1) 



The last factor being the fusion cross section is given in Section 2, and the values of x\ 2 
are found from the kinematics as 

-^quarkonium 



£1,2 



4-10- 4 e ±y . 



(1.2) 



It means that for the pp collision energy -y/s = 8 TeV, an LHCb experiment carried out 
in the rapidity range Y = 2 — 5 is in a position to measure gluon distribution with x as 
small as the record 2.5 • 10 -6 , if the lightest rj c meson is used. 



1 If one does not sum over transverse momenta of the produced quarkonia there is, strictly speaking, no 
factorization even in the LO, and Eq. (1.1) is replaced by a more complicated expression [8]. 
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The C-even quarkonia production is not the only way to probe low x partons. One can 
measure the PDFs at low x at the LHC by observing different low- mass systems, such as the 
Drell-Yan lepton pairs, or open heavy-quark QQ states. The advantage of the quarkonia 
is their direct coupling to gluons already in the LO. In the case of charmonia with their 
low mass ~ 3 GeV one achieves almost the lowest possible scale where one can justify the 
notion of the gluon distribution itself, and the use of perturbative QCD. In fact, it is not 
altogether clear beforehand if the gluon distribution at a relatively low scale corresponding 
to the charmonia production as measured in the pp collisions is not affected by power 
corrections such as the absorptive effects and/or the multiple gluon rescattering, and is 
not different from that measured, say, in the ep collisions at the same value of Bjorken x. 
This important question has to be answered experimentally. The bottomonia production 
corresponds to a higher scale, and there is most probably no such problem there. Therefore, 
comparing the gluon distributions obtained from charmonia and bottomonia production 
one would be able to judge about the possible nonlinear effects of the gluon self-interactions 
at a relatively low scale. 

There is also a theoretical problem with the simple LO Eq. (1.1). Supposing the 
inclusive b- or c-quarkonium production cross section is measured - to what precisely fac- 
torization scale hf does the gluon distribution correspond when extracted from Eq. (1.1)? 
This is an important question since one expects that the PDFs depend strongly on the 
choice of /ii? at low x because of the strong gluon bremsstrahlung there. 

In general, after summing up all orders of the perturbation theory, the final result 
should not depend on the choice of \ip that is used to separate the incoming PDFs from 
the hard matrix element a. Contributions with low virtuality, Q 2 < fj, F , of the incoming 
partons are included into the PDFs, while those with Q 2 > fi F are assigned to the matrix 
element. However, at low x the probability to emit a new gluon in an interval A/ip is 
enhanced by the large value of the longitudinal phase space, that is by the large value of 
ln(l/x). In fact, the mean number of gluons in the interval A In \lf is 



leading to the value of (n) up to about 8, for the case ln(l/x) ~ 8 and the commonly 
practiced \xf scale variation from ji/2 to 2/j,. Meanwhile, the next-to-leading (NLO), coef- 
ficient function (the hard matrix element), allows, by definition, the emission of only one 
additional parton. Therefore one cannot expect here a compensation between the contri- 
butions coming from the PDF and from the coefficient function. To that end one would 
need in this case to calculate the hard matrix element to the eighth order, which is not 
practical. [At large x the compensation is much better and provides reasonable stability 
of the predictions with respect to the variations of the scale fiF-} 

To circumvent this difficulty and to fix the factorization scale in Eq. (1.1), we use an 
approximate method following the recent Ref. [9]. The method is recalled in Section 3 
where we also find the best choice of the factorization scale hf = for the processes at 
hand: actually it determines at what scale parameter the gluon distributions are evaluated 
when the quarkonia production is measured. In the forthcoming Section 2 we evaluate the 




(1.3) 
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cross sections for the elementary hard two-gluon fusion processes into C-even charmonia. 
It becomes possible after we go into some details of the inverse processes i.e. the decays of 
charmonia. In Section 4 we discuss the resulting inclusive production cross section of el- 
even charmonia, and the ways to experimentally detect them. We stress that the absolute 
normalization of the gluon distribution obtained from the measurements we suggest, can 
be found even in the case when the experimental and/or theoretical normalization of the 
cross sections is poorly known. In Section 5 we discuss briefly the production of the 
bottomonium Xb2- The cross section is less than in the case of the charmonia production, 
however it can give an important independent information on the gluon distribution at a 
larger normalization scale. We summarize in Section 6. 



2 Elementary cross sections of the r\ c and Xco.2 production 

In the literature, one can find the LO two-gluon fusion cross sections gg — > M as well as 
the NLO differential cross sections gg —> M + g and gq — > M + q, expressed through the 
charmonia radial wave function at the origin Rq (for the s-wave charmonium r] c ) or the 
derivative at the origin R[ (for p-wave charmonia Xc0,2)- hi particular, the LO two-gluon 
fusion elementary cross sections are [10-12] 

6- LO (gg^Xco) = l2ir 2 a 2 s -^-, (2.2) 

XcO 

R' 2 

» LO (gg -> Xc2) = iQ* 2 a 2 (2.3) 

Xc2 

The numerical values of the quantities Rq and R[ have been evaluated in the past by 
many authors in the nonrelativistic quark models. Depending on the details of the model 
used these quantities lie in the ranges Rq ps (0.5 — 1.0) GeV 3 and R 2 « (0.07 — 
0.14) GeV 5 . One can try to avoid model-dependent estimates and reduce the uncertainties 
in the couplings of two gluons to the charmonia by using the experimentally-known partial 
widths of the charmonia decays. We write the C — > 77 and the C — >■ gg widths (including 
the 1st order QCD radiative corrections) from Ref. [13]: 



M„ . V \ 3 3 / vr 



r(r? c -)• 77) = ^ Qt<m ifMi+l---)^)- 



2 f 2 

T(Vc ^99) = ^ 4vr«2 (2.5) 



M Vc \ 'it 



M Xc0 V V 3 9 J vr 



2 f 2 



9 M XcQ V 



r(x c o -^99) = ^ 1 + 8 - 77 — > ( 2 - 7 ) 
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2 

Xc2 



Xc'2 



'Xc2 



L Xc2 



1 



3 7T 

4.827^ 

7T 



(2.8) 
(2.9) 



where , f Xc0 f Xc0 are the relativistic matrix elements of the local heavy-quark currents 
creating (or annihilating) the appropriate mesons from the vacuum. In the nonrelativistic 
limit they are related to the wave functions at the origin [13]: 



f 2 



3 Rl 



7T M, 



f 



108 EL? 



Vc 



XcO 



IT M 3 

XcO 



(2.10) 



with f Xc2 = f Xc0 - The system, however, is not fully nonrelativistic, and we relax this 
condition. 

We now fit the experimentally-known widths by Eqs. (2.4-2.9). We identify the gg 
widths with the total hadronic widths, T(r] c -)• gg) « r(r? c ) tot , r(x c0 -)• gg) « r(xco)tot • 
(1 — 0.0117), T(xc2 — > gg) ~ r(Xco)tot • (1 — 0.195) where in the parentheses be subtract the 
branching ratios of the radiative decays, Br(xc0 jJ/ip) = 0.0117 and Br(xc2 7J/1P) = 
0.195. Here and below the experimental numbers are from the latest PDG listings [14]. We 
treat f Vc , f XcQ , f Xc2 and a s as free fitting parameters. The results of the fit are presented 
in Table 1, and are impressively good. 





r fit (77), keV 


r exp er(77)> keV 


T m (gg), MeV 


rexper(hadrons), MeV 


Vc 


5.3 


5.3 ±0.5 


29.7 


29.7 ± 1.0 


Xo 


2.3 


2.3 ±0.23 


10.3 


10.3 ±0.6 


X2 


0.55 


0.51 ±0.043 


1.48 


1.59 ±0.11 



Table 1. A simultaneous fit to the radiative and to the hadronic widths of the C-even charmonia, 
Eqs. (2.4-2.9). 

We find the best-fit values f Vc = 432 MeV, f Xc0 = 240 MeV, f Xc2 = 361 MeV and 
a s = 0.335. Using these values in Eq. (2.10) and Eqs. (2.1-2.3) we obtain the elementary 
gluon-fusion cross sections 

a(gg ^ Vc ) ~ M4nb, (2.11) 

a(gg^ X co)-62nb, (2.12) 

a(gg^ X c2)^U0nb. (2.13) 

It should be kept in mind, though, that the QCD radiative corrections and the relativistic 
corrections to the charmonia decays appear to be rather large, therefore, the above cross 
sections extracted from the fit to the charmonia widths carry theoretical uncertainties. An 
estimate of the corrections shows that Eqs. (2.11-2.13) may be correct up to a factor of 
two in either direction. 

Indeed, the LO cross sections of the hard gluon fusion to charmonia can be derived 
alternatively from simple arguments. The decay of a spin-zero meson into two on-mass-shell 
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gluons is described by only one helicity amplitude, call it A. In terms of this amplitude 
the width of a meson with mass M is T(M — > gg) = M whereas its production cross 
section is a LO = ■ We account here for the fact that the standard gluon PDF already 
includes the sum over the 8 gluon colours and over 2 transverse polarizations. A similar 
relation between the two-gluon fusion cross section and the two-gluon decay width exists 
for the spin-2 \c2 meson; the only difference is the spin factor (2 J + 1). Therefore, the 
cross sections of the hard subprocesses can be written as 

H99 -> Vc) - ^g^^ - 539 nb ' ( 2 - 14 ) 

Vc 

o{99 -> Xco) - ^rp 124 nb, (2.15) 

XcO 

5n 2 T(Xc2 99) 

o{gg -> Xc2) ^3 ~ 85 nb , (2.16) 

where for the numerical evaluation we have replaced the two-gluon widths by the phe- 
nomenological hadronic widths as above. 

Comparing the estimates (2.11-2.13) with the estimates (2.14-2.16) one gets the idea of 
the theoretical uncertainty in evaluating the elementary cross sections. The first derivation 
takes into account the radiative corrections to the charmonia decays but ignores them in 
the cross sections. The second derivation is based on the fact that the effective Cgg vertex 
(C = Xc Vc) is the same in the decay into two on-mass-shell gluons as in the fusion of two 
gluons (that should be considered as being on-mass-shell in the LO) into a charmonium. 
In both cases the radiative corrections seem to be the same. Therefore, we are inclined to 
trust more the second estimate (2.14-2.16), given the experience of the first one: it shows 
that the two-gluon decays can be well replaced by the total hadronic widths. 

3 The scale parameter for the gluon distribution 

To sketch the idea how to choose the appropriate scale, we start with the LO expression 
for the cross section. In the collinear approach, the cross section has the form 

a{fx F ) = PDF( / u F )®C L °^PDF(/i F ), (3.1) 

where C LO denotes the LO hard matrix element squared. The effect of varying the scale 
from m to //p in both PDFs can be expressed, to the first order in a s , as 

<r(ji F ) = PDF(m) ® (c LO + gin (&\ (Pi cft C LO + C LO P rig ht)) ® PDF(m), (3.2) 

where the splitting functions Pi c f t and P r i g ht act on the left and on the right PDFs, respec- 
tively. Let us recall that in calculating the a s correction in Eq.(3.2), the integral over the 
transverse momentum (virtuality) of the parton in the LO DGLAP evolution is approxi- 
mated by the pure logarithmic dk 2 /k 2 form. That is to say, in the collinear approach, the 
Leading Log Approximation (LLA) is used. 
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Let us now study the cross section at the NLO. First, we note that the original Feynman 
diagrams corresponding to the NLO matrix element C N formally do not depend on \ip. 
However, we shall see below that in fact scale dependence appears. In the NLO we can 
write 

a(n F ) = PDF(hf)® (C LO + a s C c ^ r °)«)PDF(/i F ), (3.3) 

where we include the NLO correction to the coefficient function. In terms of Feynman 
diagrams it means that the gg — > i] c (Xc) subprocess plus the 2 — > 2 subprocesses, gg — > 
r lc(Xc) + 9 and qg — > r] c (xc) + 9> ar e now calculated with better than the LLA accuracy. 
However part of this contribution is already included, to the LLA accuracy, into the second 
term in Eq.(3.2). Therefore this part should be now subtracted from C NL °. Moreover, this 
LLA part depends on the scale hf- As a result, changing hf redistributes the order a s 
correction between the LO part (PDF ® C LO ® PDF) and the NLO part (PDF <g> a s C^? <g> 
PDF). 

We see that the part of the NLO correction that remains after the subtraction, C^P(hf), 
depends now on the scale hf as due to the hf dependence of the LO LLA term that has 
been subtracted out. The trick is to choose an appropriate scale hf = ho such as to 
minimize the remaining NLO contribution C^ m (^ F ). To be more precise, we choose the 
value hf = Ho such that as much as possible of the 'real' NLO contribution (which has a 
ladder-like form and which is strongly enhanced by the large value of ln(l/x)) is included 
into the LO part where all the logarithmically enhanced a s ln(l/x) terms are naturally 
collected by the incoming parton distributions 2 . 

As shown in Ref. [9], after the scale hf = Ho 1S fixed for the LO contribution the varia- 
tion of the scale in the remaining NLO part does not change noticeably the predicted cross 
section. Moreover, it was shown that in the case of the Drell-Yan lepton pair production 
the NLO prediction with hf = Ho is very close to the NNLO result. 

We now determine the "best" value of the scale hf = Ho f° r which the factorization Eq. 
(1.1) is maximally correct. It will be in fact the scale parameter for the gluon distribution 
measured from the rj c , Xc inclusive production rates, if one uses Eq. (1.1) to determine the 
gluon distribution. 

As explained above, in order to find the value of the appropriate scale hf = Ho °f the 
LO contribution we have to know the cross section of hard subprocesses calculated at the 
NLO level. The differential cross sections of the gg — >■ M + g and gq — >• M + q subprocesses 
as functions of the Mandelstam variables s and t are presented in Refs. [12, 17, 18] and 
are collected in the Appendix. We integrate them there over the available t interval, 
subtract the contributions generated by the last step of the LO DGLAP evolution up to 
the factorization scale hf, convoluted with the LO cross sections. Finally, we choose the 
value of the factorization scale hf — Ho such that it nullifies the remaining NLO gg —> M+g 
and gq — ^ M + q contributions. 

2 Actually our approach is rather close in spirit to the fct-factorization method. Using the known NLO 
result we account for the exact kt integration in the last cell adjacent to the LO hard matrix element 
(describing the gg — > 7? c (Xc) boson fusion), while the unintegrated parton distribution is generated by the 
last step of the DGLAP evolution, similarly to the prescription proposed in Refs. [15, 16]. 
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The value of found by this method may be in fact different for various subprocesses. 
It depends also on the subprocess energy s. Therefore we have to average the gg — >■ M + g 
and gq — >■ M + q cross sections with the incoming parton flux F(s) driven by the PDF low 
x behaviour. For low- a; parton distributions, we assume a power behaviour, F{s) oc s _A 
with the power < A < 0.3. Depending on the choice of A, we present in Table 2 the scale 
V^F = ^0 * na ^ nullifies the remaining NLO contribution of the gq — > rj c + q, gq — > Xc(J) + 1 
and of the gg — > rj c + g, gg — > Xc(J) + 9 subprocesses. 



subprocess 


A = (s ->■ oo) 


A = 0.1 


A = 0.2 


A = 0.3 


gq^-rjc + q 


3.3 


3.1 


2.9 


2.75 


gq -> x(o) + 9 


2.4 


2.3 


2.2 


2.1 


59 -> X(2) + 9 


2.9 


2.8 


2.6 


2.5 


99 ^-Vc + 9 


3.3 


3.0 


2.75 


2.5 


99 -> X(0) + 9 


2.4 


2.1 


1.9 


1.7 


99 -> X(2) +9 


2.9 


2.5 


2.2 


1.9 



Table 2. The best scale /Ug (in GcV 2 ) calculated from various subprocesses, depending on the 
power A in the gluon flux assumed. 

In the case of an asymptotically high subenergy s — > oo, when the ladder-type diagrams 
dominate, the values of /io are the same for both gq — >■ M + q and gg — > M+g subrpocesses. 
However even the LHC energy is not sufficient to reach the asymptotics. Actually the 
rapidity interval available at the LHC, SY ~ 10, corresponds approximately to A ~ 0.1. 
This value looks also as realistic for the gluon distribution at low x and relatively low scale 
~ 2.5 GeV 2 . For A > the value of fio needed to nullify the remaining NLO contribution of 
gq —> M + q subprocess is larger than that for the gg —> M + g case. Let us note, however, 
that in the last case by changing the value of \ip we try to mimic by the LO-generated 
contribution also the terms that have the structure rather different from that generated by 
the LO evolution. This is not altogether consistent. Therefore we believe that the value 
of fiQ calculated from the gq —> M + q subprocess whose Feynman diagram has the same 
form as that generated by the DGLAP evolution, is more reliable. 

It is interesting that for the pseudoscalar rj c production we get a larger value of no 
despite that its mass is less than that of Xc0,2- Owing to the unnatural parity of r] c , 
the production vertex contains an additional transverse momentum that enhances large- \t\ 
contributions. To compensate it, one has to take a larger fiQ. 

We see from Table 2 that we still have some ~ 10 — 20% uncertainty in the value of 
the appropriate scale /io but this is much less than the usually used ad hoc interval from 
M/2 up to 2M. Moreover, when and if it comes to fitting the data it will be possible to 
simultaneously specify/determine the value of A and to fix the appropriate scale (ip = W) 
more precisely. At the moment we think that the power A = 0.1 is the most realistic for a 
relatively low scale fi F = 2 — 3 GeV 2 . 
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4 Extracting gluon distribution from the C-even charmonia production 



Choosing the appropriate scale \ip = (J,q from Table 2 we strongly suppress the remaining 
higher a s order contributions to the LO factorization Eq. (1.1). Thus, the inclusive cross 
section of the r] c and Xc0,2 production 



da 
dY 



x±g(xi, no) X2g(x2, Ho) &{gg — > C— even charmonium) 



(4.1) 



will measure directly the product of gluon densities at the normalization point fio- The 
values of x\ % are found from Eq. (1.2) while the values of a are given by Eqs. (2.14-2.16). 

For example, at yfs = 8 TeV and Y = 5 we have for the r\ c production x\ = 0.055, X2 = 
2.5 • 10 -6 , and for the Xc2 production x\ = 0.066, X2 = 3.0 • 10 -6 . At the upper side (xi) 
the gluon distribution is rather accurately established, therefore from measuring the rate 
of the charmonia production one can find the gluon densities at unprecedented low X2, see 
in more detail below. 

Let us briefly discuss how to register the production of the r\ c and Xc mesons. Of the 
three C-even charmonia considered the most favourable observational conditions seem to be 
for the Xc2 meson via an anomalously large radiative decay Bv(xc2 lJ/ip) = 0.195±0.008. 
Actually the Xc2 inclusive production has been already observed at the LHCb via this 
particular decay channel [19]. The expected Xc2 production rate, times this branching 
ratio, times the branching ratio Br(J/ip — > [i + = (5.93 ±0.06) • 10~ 2 is plotted in Fig. 2. 
It appears to be quite large - in the range of tens of nanobarns. 

do-/dY * Br * Br, nb 

30- 



Figure 2. Cross section of the inclusive \ci production per unit rapidity Y, times the branching 
ratio of its decay into jj/tp, times the branching ratio of the J/ip — > [i + H~ decay, in nanobarns. 
The shaded area corresponds to the gluon flux from the shaded area in Fig. 1, whereas the dot- 
dashed, solid and dashed curves correspond to the extrapolation using x g(x) ~ const. (A = 0), 
xg(x) ~ 1/x 01 (A = 0.1), and xg(x) ~ l/x - 2 (A = 0.2), respectively, shown in Fig. 1, left. The 
scale parameter /ig = 2.5 GeV 2 is assumed for the gluon distribution. 

The XcO meson has a comparable production rate but a much smaller radiative decay 
branching ratio Br(xc0 ~~ ^ lJ/ip) = 0.0117± 0.008. The production times branching curves 
for the XcO meson are similar to those shown in Fig. 2 but the overall scale is an order of 
magnitude less. Therefore, unless a good hadronic decay channel is found, the XcO meson 
cannot compete with its Xc2 cousin. 
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Finally, the r\ c meson decays mainly into it, K mesons and it is not easy to find a meson 
decay mode that would not be plagued by the huge multi-meson combinatorial background 
at the LHC energies. Probably the best channel would be the decay into pp with the 
branching ratio Br(?7 c — > pp) = (1.41 ± 0.17) • 10~ 3 or into AA with the branching ratio 
Br(?7 c AA) = (0.94 ± 0.32) • 10" 3 . The r] c cross section times the pp branching is similar 
and close in magnitude to what is presented in Fig. 2 for the \ci meson. It should be 
noted that the r\ c production measures, paradoxically, the gluon distribution at a larger 
scale than the Xc2 meson (see Section 3) and therefore the two measurements can be of 
independent value. 

It is interesting that when in the whole x region probed by an experiment the gluon 
distribution has the power-law form g(x) ~ l/x 1+A the production cross section is inde- 
pendent on the rapidity Y of the charmonia produced. However the hight of the plateau 
is extremely sensitive to the power A at small x. 

Remarkably, even if one does not know the absolute normalization of the experimental 
cross section for the charmonia production and/or of the theoretical cross section (4.1), one 
can still extract the absolute normalization of the gluon distribution at low x, by matching 
the measurements with the distribution in the medium- x range where it is already known. 

To give the idea, let us consider the setup of the LHCb where particles with rapidity 
up to Y = 5 can be registered. Supposing it is found that the number of counts C(Y) 
of the Xc2 mesons produced in a rapidity interval Y ± 6Y is roughly independent of Y, 
corresponding to the plateau in Fig. 2. It means that the gluon distribution at very small 
x has a power-law behaviour, xg{x) = ax~^. We want to find the power A and the 
absolute normalization a. We write 

x 1 g(x 1 )x 2 g(x 2 ) = MC{Y), x li2 = ^e ±y , (4.2) 

V s 

2/_*aj = tfC(Y), (4.3) 

where J\f is an unknown normalization factor 3 . We see from Eq. (4.3) that the power law 
for g(x) is the only one leading to the number of counts C(Y) independent of the rapidity. 

However, at Y > 4 the number of counts will deviate from the plateau, see Fig. 2. 
This is where one of the two fusing gluons has a known flux xig(x±), see the solid part 
of the curve at the right-hand side of Fig. 3. At the end point of the solid curve, that is 
at x = x c ~ 0.028 the gluon distribution is still known to an accuracy of a few percent . 
Therefore, one can determine the other gluon's xig{xi) from Eq. (4.2): 

1 , A C ( Y ) a -A ,x 

A7 X2 ^ 2) = ^R)=A7^ • (44) 

To find numerically A and fy, one has to make two or more measurements at Y > 4, for 
example as shown in Fig. 3, and make a two-parameter fit to Eq. (4.4). Then, assuming 



3 In fact TV -1 is the elementary fusion cross section (2.13), times the integrated luminosity, times the 
registration efficiency. 

4 The accuracy can be estimated by comparing two known gluon distributions [2-4] at x c . 
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x g(x) 

1 8 r 

Y=4.25 




Figure 3. Finding the gluon distribution xg(x) from Eq. (4.2). The solid part of the curve shows 
the known distribution at the normalization point fj,Q = 2.5 GeV 2 . The dashed part is the supposed 
power law at very small x, xg(x) = ax~ A . The rectangular "gates" indicate the values of x\ 
(right point) and X2 (left point) in the product of the gluon distributions, corresponding to a given 
rapidity Y of the y c2 meson produced, y/s = 8 TeV is assumed. 

the power behaviour of the gluon distribution all the way up to x c we equate 

jfX~ A = ^x c g(x c ). (4.5) 

In this equation, the combination jj- is presumably known from the fit above, and all the 
rest quantities are also known, except the normalization factor N . Therefore, Eq. (4.5) 
enables one to find M and hence the absolute normalization of the gluon distribution a. 

Alternatively, one can find A, a and J\f separately by solving the system of equation 
(4.3) and two equations (4.4) evaluated at two different rapidities Y > 4. 

If the actual behaviour of the gluon distribution at very low x is substantially different 
from power-like, this will be seen from the deviation from the flat plateau in the production 
rate as function of Y. The data should be then analyzed accordingly, however in any case 
the absolute normalization of the gluon distribution will be possible to deduce even without 
knowing the absolute values of the charmonia production cross section - by matching the 
data with the gluon distribution at x > x c ~ 0.028 where it is already known with a 
reasonable accuracy. 

We would like to remark that a good complement would be measuring charmonia 
production in a fixed-target experiment with the LHC beams (AFTER@LHC) as it allows 
to observe charmonia with low pr and to extract the gluon distribution at x from a few 
units of 10" 3 to x ~ 1 [20, 21] 5 . 

5 Gluon distribution from the C-even bottomonia production 

The same theoretical considerations can be applied to measuring gluon distributions from 
the production of the C-even bb mesons, such as Xfc2(l-P)(2 ++ , 9912). Like Xc2, the bot- 
tomonium Xb2 has a large branching ratio for the radiative decay, Br(x62(l-P) — > 7^(15)) = 
0.191 ± 0.012 while the leptonic branching ratio for the T is Br(T(15) — > fi + fi~) = 

5 We thank J. -P. Lansberg for bringing our attention to this work. 
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(2.48 ±0.05) • 10~ 2 . This decay cascade makes the observation of the bottomonium Xb2^-P) 
possible. 

At large heavy-quark masses, the gluon-fusion cross sections of the quarkonia produc- 
tion scale as cr(gg — > x) ~ a 1{^x) I anc ^ therefore the Xb cross section is expected to be 
~ 20 times less than that of Xc- Prom the evaluation of the Xc2 production cross section 
(2.16) we estimate a(gg — > Xb2i)-P)) ~ 4nb. The Xb2(2P) production cross section must 
be 3 — 4 times smaller, according to the nonrelativistic estimate of R'i(0). 

However, the smaller elementary Xb2 production cross section is multiplied in Eq. (4.1) 
by a larger gluon flux x\ g(x±, /Uq) X2 g(x2, Mo) expected at the scale /1q appropriate for the 
bottomonia as contrasted to the charmonia. According to the derivation in Section 3, the 
scale //q is proportional to the mass squared of the quarkonium in question. It is known 
that at higher resolution scale the gluon distribution increases towards small x effectively 
as a higher power x~ A . Taking A ~ 0.25 and using Table 2 we find that for the Xb2 
production the gluon distribution scale is /x 2 , ~ 20GeV 2 ; it is plotted in Fig. 1, right. 

dtr/dY * Br * Br, nb 

2.0 r 




Figure 4. The expected cross section of the inclusive bottomonium Xb2(lP) production per unit 
rapidity Y, times the branching ratio of its decay into jT(1S), times the branching ratio of the 
T(liS') — > /i + fi~ decay, in nanobarns. The gluon distribution shown in Fig. 1, right, is assumed. 

The expected cross section of the inclusive bottomonium XbzQP ) production, times the 
branching ratio of its decay into "fT(lS), times the branching ratio of the T(1S") — > 
decay, is plotted in Fig. 4, right. It is not a plateau anymore, even if one assumes a power- 
law behaviour of the gluon distribution at low x. We see that it is a few times less than the 
inclusive Xc2 production, times its branching ratios, see Fig. 2, but probably within reach. 

The bottomonium Xb2(2P)(2 ++ , 10269) can be observed via the radiative decay into the 
twoT's with the branching ratios Br(x& 2 (2P) ->■ jT(lS)) = 0.071±0.01 and Br(xfe 2 (2-P) ->■ 
7T(25)) = 0.162 ± 0.024. The consequent leptonic branching ratios for the T decays are 
Br(T(15) -> /U+aO = (2.48 ± 0.05) • 10~ 2 and Br(T(25) -> /i + /x - ) = (1.93 ± 0.17) • 10~ 2 , 
respectively. Combining these decay cascades, the total registration rate of the Xb2(^P) 
bottomonium is expected to be very similar to that shown in Fig. 4, right, being however 
3 — 4 times less. 

Although the production cross section (times the branching ratios) for Xb2{^P) is 
several times less than that for Xc2(^-P), it would be easier to match the measured gluon 
distribution at very low x to that already known at larger values of x, see the end of Section 
5. The smallest value of x accessible from the bottomonium production is x m ; n = 8.3-10" b . 
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It should be noted that the Xb2 must have a broad distribution in the transverse momenta 
as due to the typical double- logarithmic QCD form factor [8]. 

6 Conclusions 

The inclusive production cross sections of C-even charmonia Xc2{2 ++ , 3556) and r] c (0 \ 2998) 
at the LHC, times the branching ratios of their convenient decay modes, Br(x C 2 — > iJ/i 1 ) • 
Br( J/ip — > ^ + /i~) and Br(7y c — > pp), respectively, are estimated to lie in the range 5 — 30 nb, 
depending on what is the actual behaviour of the gluon distribution at very low Bjorken x. 
Measuring the production of those charmonia, integrated over their transverse momenta 
will enable one to determine the fundamental quantity — the gluon distribution in nucleons 
g(x,Q 2 ) — at an unprecedented low x > 2.5 • 10 -6 and relatively low normalization scale 
Q 2 = 2.5 — 3Gev 2 . The absolute normalization of the gluon distribution can be found 
by matching the measured charmonia yield with the gluon distribution at higher x where 
it is already known, even if the normalization of the experimental and/or theoretical cross 
sections are not well established. 

Similarly, measuring the inclusive production of the bottomonium Xb2(2 ++ , 9912) with 
the cross section times the branching ratios around 1 nb will allow to extract the gluon 
distribution at x > 8.3 • 1CT 6 but a larger scale Q 2 « 20GeV 2 . 

Combining the measurements of the two quarkonia production will give a rather full 
knowledge of the fundamental quantity - the gluon distribution - in a broad range of x 
and Q 2 . In particular, one will be in a position to judge if at Q 2 = 2.5 — 3Gev 2 the 
nonlinearity (the gluon self- interaction) becomes important or not, and to discriminate 
between various theoretical models of the high-energy processes. 
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A Appendix 

We list below the LO and NLO cross sections of r] c and Xc mesons production [12, 17, 18]. 
We have checked the equations and present them in the form which is further on used for 
integrating the differential cross sections over t. 

The LO contributions to the hard gg — > Xc(J = ®, 2) and gg — > r] c (J = 0) subprocesses, 
i.e. the LO two-gluon fusion cross sections are given by Eqs. (2.1-2.3). 

The NLO 2—7-2 differential cross sections are expressed through the quantity s denot- 
ing the x + X + Q or 7 ? + 5) V + Q energy squared, r = s/M 2 . The variable z is defined by 
the equation t = — s(l — 1/r + z)/2, where t is the momentum transfer squared from the 
initial to the final gluon or quark. 
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The NLO cross section is obtained by integrating the differential cross section over 
t, that is translated in our notations into the integration over z. The divergencies of the 
integrands at z = 1 — 1/r or at z = — 1 + 1/r reflect the logarithmic divergencies of the 
differential cross sections at t — > or at u = M 2 — s — t — > 0. They are the collinear 
singularities that are responsible for the evolution of the PDF, 

^ I = -La LO K(l/r)r 

where K(x) is LO DGLAP splitting function for the gg or the qg channels. To avoid 
the double counting we subtract from the NLO cross sections the logarithmic part at 
z > 1 — 1/r — 2fj, F /s (and at z < — 1 + 1/r + 2fj, F /s if there is a singularity in u) as being 
attributed to the PDF, thereby removing the infrared divergency at t — > 0. 

Since the NLO qg — > qM cross section is described by the same diagram as that 
responsible for the LO DGLAP evolution we choose the scale \xf = A*o such that being 
integrated up to /io the LO-generated contribution nullifies the remaining NLO qg — > qM 
cross section. By doing that we shift the major part of the corrections (enhanced by the 
large value of ln(l/x)) to the low x parton distributions. Below we list the NLO cross 
sections used in this derivation. 

The gg — > g + Xc(0) differential cross section is 

^^9 + *M) _«^ tueu F x (0) = Ar x (0)/D x s(0) (A.l) 

Xc(0) 

where 

AT* (o) = -32 [3(154 r + 27) + {z 2 + 3) 4 (z 2 - 1)V 4 - 12(76z 2 + 159)r 3 

-(270z 2 - 187)r 2 + 2(87z 4 + 848z 2 + 649)r 5 + (279z 4 + 1004z 2 - 663)r 4 
-2(3z 6 + 25z 4 + 85z 2 + 47)(z 2 + 3) 2 (z 2 - l)r 13 

-(36z 6 + 893z 4 + 3418z 2 - 4875)r 6 + 8(42z 6 + 517z 4 + 1268z 2 - 507)r 7 

+2(69z 6 - 531z 4 - 6217z 2 - 5193)(z 2 - l)r 9 

-(81z 8 + 1760z 6 + 5858z 4 + 16312z 2 + 7669)r 8 

+(18z 10 + 12492 8 + 10424/ + 34958z 4 + 21726z 2 + 2025)r 10 

-4(48z 10 + 661z 8 + 3172z 6 + 7958z 4 + 5604z 2 + 1757)r n 

+(9z 12 + 212^0 + 1809z 8 + 5952z 6 + 11019z 4 + 5516z 2 + 1083)r 12 ] , 

D* g (0) = {rz + r + l) A {rz + r - l)(rz - r + l)(rz - r - l) 4 (r - l) 4 r . 
We integrate it over t and obtain 

- NL %9^9 + X c m = ^^^ g (0) (A.2) 

Xc(0) 
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where 



s — M 2 

2* (0) = [64((1081n ^(r 2 - r + l) 2 (r + l) 4 (r - l) 2 

-(172r 10 - 56r 9 - 617r 8 + 188r 7 + 1104r 6 - 508r 5 + 302r 4 
+52r 3 + 92r 2 + 132r + 99))(r 2 - 1) 

-12(9r n - 31r 9 + 14r 8 + 40r 7 - 10r 6 - 176r 5 + 42r 4 + 7r 3 
+10r 2 - 41r - 24) ln(r)r)] / [3(r + l) 5 (r - l) 4 r 2 ] . 

The <7<7 g + Xc(2) differential cross section is 

gto^^g)) _ JpgL^ (2) , F * (2) = „* (2)/£) * (2) (A . 3) 

Xc(2) 



where 

iV* (2) = 64[6(28r + 9) + (z 2 + 3) 4 (z 2 - l) 2 r 14 

+6(5z 2 - 34)r 3 - 5(36z 2 + 103)r 2 + 2(93z 4 + 1405z 2 + 2022)r 4 

-2(519z 4 + 1567z 2 + 2762)r 5 - 2(3z 6 - 32z 4 + 199z 2 - 10)(z 2 + 3) 2 (z 2 - l)r 13 

-(24z 6 + 3017z 4 + 10102z 2 + 3705)r 6 + 4(315z 6 + 2264z 4 + 5635z 2 + 2154)r 7 

-2(27z 8 + 178z 6 + 760z 4 + 7742z 2 + 605)r 8 

-4(93z 8 + 1197/ + 4055z 4 - 281z 2 + 216)r 9 

+ (12z 10 + 1099z 8 + 8732z 6 + 29186z 4 + 6336z 2 - 309)r 10 

-2(21z 10 + 482z 8 + 4190z 6 + 13432z 4 + 837z 2 + 1006)r n 

+2(3z 12 - llz 10 + 186z 8 + 2634z 6 + 5991z 4 - 2927z 2 + 780)r 12 l , 



D* (2) = 3(rz + r + l) 4 (rz + r - l)(rz - r + l)(rz - r - l) 4 (r - l) 4 r . 
We integrate it over t and obtain 



where 



s — M 2 

2* (2) = [I28((721n ^(r 2 - r + l) 2 (r + l) 4 (r - l) 2 

t^F 

-(106r 10 - 32r 9 - 101r 8 + 239r 7 + 651r 6 - 793r 5 + 395r 4 
-527r 3 + 35r 2 + 201r + 66))(r 2 - 1) 
-12(6r n - 22r 9 + 8r 8 - 74r 7 - 31r 6 

-llr 5 + 204r 4 - 86r 3 - 17r 2 - 5r - 12) ln(r)r)] / [9(r + l) 5 (r - l) 4 r 2 ] . 

The gg —> g + i] c differential cross section is 

dajgg -» g + V ) vra|R| 

dt 4s 2 M| 99 991 99 { 1 
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where 

N 71 = 

99 



-{r 4 z 4 + 6r 4 z 2 + 9r 4 - 12r 3 z 2 

-4r 3 + 6r 2 z 2 + 6r 2 - 4r + 9)(rV + 3r 2 - 2r - l) 2 , 



"^33 = (rz + r + l) 2 (rz + r — l)(rz — r + l)(rz — r — l) 2 (r — l) 2 r . 



We integrate it over t and obtain 



8sM 3 99 



where 



s — M 2 

Tl= [ 2 ((121n^^(r 2 -r + l) 2 (r + l) 2 

-(12r 6 + 23r 4 + 24r 3 + 2r 2 + ll))(r 2 - 1) 

-12(r 7 - 5r 5 - 2r 4 - r 3 - 3r - 2) ln(r)r)] / [3(r + l) 3 (r - l) 2 r 2 ] 
The qg ^ q + Xc(0) differential cross section is 

"< m ^ + »-<°» = l^^(O), F* (0) . iV* (0)/ D * (0) 

Xc(0) 



where 



99 

We integrate it over i and obtain 



iV* (0) = (r 2 (z 2 - 2z + 5) + 2rz -2r + l)(rz + r + 5)' 
D*J0) = 2(rz + r+ l) 4 (rz + r - 1) . 



Xc(0) 



where 



s-M 2 

T* (0) = [27 Qn ^(2r 2 - 2r + l)r 

-2(43r 2 - 14r + 4)(r - 1) - 6(9r 2 - 9r + 4) ln(r)r] /(6r 2 ) . 
The qg — > g + Xc(2) differential cross section is 

?'2 



g<*^|±*^> - W W = (2)/ D * (2) 

where 



iV* (2) = [r 4 4 + 2z 2 + 8z + 5) - r 3 (44z 2 + 8z - 36) + r 2 (22z 2 - 48z + 

+48rz - 4r + 25] , 
£>* (2) = (3(rz + r + l) 4 (rz + r - 1)) . 



We integrate it over t and obtain 



* NL ° ( q9 ? + Xc(2)) = ^S^^(2) (A.10) 



where 



T* (2) = [18 In ^(2r 2 - 2r + l)r 

-(53r 2 - 16r + 20) (r - 1) - 3(12r 2 - 12r + 5) ln(r)r] /(9r 2 ) . 
The qg ^ q + f] c differential cross section is 

* (w " + - 4TOlflS -C «, = Jva/flt (A.n; 



where 



(it 9s 2 M3 99 -191-19 



N% g = 2(z - l)r + 1 + (z 2 - 2z + 5)r 2 
m = (rz + r + l) 2 (rz + r - 1) . 



95 

We integrate it over t and obtain 

^°to^, + ,) = ^^ (A.12) 

where 

r s — M 2 „ 

T ?s= [ ln 2-^(2r 2 -2r + l)-2(ln(r) + l)(r-l)r]/r. (A.13) 

After subtracting the logarithmically divergent parts of the NLO cross sections (at- 
tributed to the PDF) we have to average the remaining cross sections over the incoming 
subenergy s, that is to integrate over r with the weight driven by the parton flux F(s). 
Assuming the power behaviour x g(x) oc of the low-x gluon distribution we obtain 

the flux F oc r~ A . Therefore, our goal is to choose such a scale /j, f = that nullifies the 
integral 

f°° dr 

*fAV^K A y. (A.14) 



The resulting scales /j,q depending on the subprocess and on the power A are presented in 
Table 2. 
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